Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 October 2010 (MN KTeX style file v2.2) 



Panchromatic Averaged Stellar Populations: Pa asp 



o 

(N 

» 

o 
O 



O 

u 

Oh 
6 

o3 



> 



o 
o 



X 



R. Riffel 1 *, C. Bonatto 1 , R. Cid Fernandes 2 , M. G. Pastoriza 1 and E. Balbinot 1 

1 Departamento de Astronomia, Universidade Federal do Rio Grande do Sul. Av. Bento Goncalves 9500, Porto Alegre, RS, Brazil- 
2 Departamento de Ftsica - CFM - Universidade Federal de Santa Catarina, Floriandpolis, SC. Brazil 



ABSTRACT 

We study how the spectral fitting of galaxies, in terms of light fractions derived in one 
spectral region translates into another region, by using results from evolutionary synthe- 
sis models. In particular, we examine propagation dependencies on Evolutionary Popu- 
lation Synthesis (EPS, grasil, galev, Maraston and galaxev) models, age, metallic- 
ity, and stellar evolution tracks over the near-UV — near infrared (NUV — NIR, 3500A to 
2.5/im) spectral region. Our main results are: as expected, young (t < 400 Myr) stellar popu- 
lation fractions derived in the optical cannot be directly compared to those derived in the NIR, 
and vice versa. In contrast, intermediate to old age (t > 500 Myr) fractions are similar over 
the whole spectral region studied. The metallicity has a negligible effect on the propagation of 
the stellar population fractions derived from NUV — NIR. The same applies to the different 
EPS models, but restricted to the range between 3800 A and 9000 A. However, a discrepancy 
between GALEV/Maraston and GRASIL/GALAXEV models occurs in the NIR. Also, the initial 
mass function (IMF) is not important for the synthesis propagation. Compared to STARLIGHT 
synthesis results, our propagation predictions agree at ~95% confidence level in the optical, 
and ~85% in the NIR. In summary, spectral fitting performed in a restricted spectral range 
should not be directly propagated from the NIR to the UV/Optical, or vice versa. We provide 
equations and an on-line form (Panchromatic Averaged Stellar Population - PaASP) to be 
used for this purpose. 
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1 INTRODUCTION 

A key issue in modern astrophysics is to understand how galax- 
ies form and evolve, and the study of the stellar population and 
star formation history (over a wide wavelength range) may provide 
clues to the dominant mechanism. For example, two main scenarios 
are proposed to explain the formation of galaxies, one is the tidal 
torque theory , which suggests an initial collapse of gas at very hig h 
redshift (e.g. lEggen. Lvnden-Bell & Sandagdfl962l : IWhitelll984l) . 
and another is ass ociated with galaxy mergers at intermediate to 
high r edshifts (see ISearle & Zinnll978l : lHarnmer et alj2005Ll2007l 
I2009L for example). Both produce distinct signatures in the result- 
ing stellar populations. 

As we still cannot access the light of individual stars in al- 
most any galaxy beyond the local group, the way widely used 
is to study the stellar populations by the integrated light of the 
whole galaxy or a significant fraction of it, by means of long-slit 
spectroscopy. The usual approaches for studying unresolved stellar 
populations are either by means of empirical popul a tion synthesis 
i Bica & Alloinlll987l; |Bicalll988t iBica et al.lfl99ll: jBonatto et al.l 



2003 |: lMarastonlll998l.l2005l ; IVazdekis & Arimotoll999l : ISilva et al.l 



19981 : iKotulla et alj|2009f). as well as a combination of both tech- 
nique s dCid Fernandes et ai]|2004 120051 : iRiffel et alj[200l 120081 
120091) . By using them, we can infer the age and metallicity distri- 
butions of the stellar populations that make up a galaxy's spectral 
energy distribution (SED). 

The study of unresolved stellar populations in galaxies, 
ranging from star-forming to ellipticals, as well as in compos- 
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other hand, addressing unresolved stellar populations in the 
near infrared (NIR ) is less common, starting ~3 decades ago 
dRieke et"al 198Ct e.g.), but which seems t o be in increasing 
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Very recently. Ichen et alj feOlCh have shown that using differ- 
ent EPS models in the optical leads to different stellar population 
results. So, it is not reasonable to directly compare stellar popula- 
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tions estimated from different EPS models. They suggest that to get 
reliable results, one should use the same EPS models to compare 
different samples. As related issues, can one compare the light- 
fraction results derived, with the same base of elements, in one 
spectral region to another? Are the population fractions derived in 
the optical the same in the NUV/NIR? What is the effect of the nor- 
malisation point used in the synthesis? Does the choice of elements 
that compose the base (i.e. age, metallicity and initial mass function 
- IMF) produce different results along the whole spectrum? 

The answer to the above questions would be easy if stellar 
population synthesis were based on a more physical parameter, 
like mass-fractions. However, given the non-constant stellar mass- 
to-light (M/L) ratio, mass-fractions have a much less direct rela- 
tion with the observables than light-fractions. In this context, we 
feel motivated to carry out a detailed study of the panchromatic av- 
eraged stellar population (PaASP) components over the 3500A to 
2.5/im spectral region, commonly used for stellar population syn- 
thesis. It should be clear that PaASP is not a tool for performing 
stellar population synthesis. Instead, it is specifically designed 
for translating the spectral fitting of galaxies (in terms of light 
fractions) derived in one spectral domain into another. This pa- 
per is structured as follows: The EPS models used are described in 
Sect. [2] The methodology and results are presented in Sect. [3] Re- 
sults are discussed in Sect.[4] The final remarks are given in Sect.[5] 



clusters to galaxies, both in ter ms of resolved stel lar populations 
and integrated light properties dKotulla et al1l2009l) . According to 
the authors, the code considers both the chemical evolution of the 
gas and the spectral evolution of the stellar component, allowing 
for a chemically consistent approach. Thus, some SSPs provided 
by GALEV show a emission line spectra. 

The SSP models provided by GALEV cover 5 metallicities and 
4000 ages (0.02 < Z/Z Q 2.5 and 4Myr < t 16Gyr). They 



are based on the spectra from BaS eL spectral library l lLeieune et al 



1 199711 19991 ; IWestera et ai]l2002l) . originally based on the 



Kurucz 



d 19921) library . They also hav e 3 different IMFs JSalpetej|l95l - 



lKroupdl200~j| ; IChabriedliool) and use the theoretical isochrones 
from the Padova team (e.g. lBertelli et alJl994l : ISchulz et al]2002ri . 
The wavelength coverage spans the range from 90 A to 160/^m, 
with a spectral resolution of 20A in the NUV-optical, and 50-100A 
in the NIR. They also include the thermally- pulsating asymp- 
totic giant branch (TP-AGB) phase pro vided by the Padova tracks 
dSchulz et ai]|2002l ; iGirardi et alj|2002h . It is also worth mention- 
ing at this point that the TP-AGB stars account for 25 to 40% of the 
bolometric light of an SSP, and for 40 to 60% of the lig ht emitted 
in the K-band (e.g. ISchulz et alj|2002l; lMarastoij|2005L and refer- 
ences therein). In addition, contrary to lMara ston (2005, see below), 
GALEV models do not include e mpirical TP-AGB sp ectra. For a de- 
tailed description of GALEV see iKotulla etal1 ( l2009t) . 



2 THE ADOPTED EPS MODELS 

In this section we describe the EPS models used in what follows. 
The four models: GRASIL, GALEV, Maraston and GALAXEV, have 
been selected because they have a spectral coverage from the NUV 
to the NIR region (~3500A - 2.5/im) and are widely used. A brief 
description of them is made below. D etails can be found in the orig- 
inal papers cited below, as well as in IChenetaljJioiol) . 



2.1 GRASIL 

The GRAphite a nd SILicate - GRASIl|3 - code, developed by 
ISilva et alj (11998), is a chemical evolution code that follows the 
star formation rate, metallicity and gas fraction, which are ba- 
sic ingredients for stellar population synthesis. The latter is per- 
formed with a grid of integrated spectra of simple stellar popula- 
tions (SSPs) of different ages and metallicities, in which the ef- 
fects of dusty envelopes around asymptotic giant branch (AGB) 
stars are in cluded, but not the AGB energetics dBressan et al.1 
1 19981 l2002h. The models consider four initial mass functions 
(IMFs). iKuruczl dl992h; iKennicutt et al.1 dl994h : IScalol < ll986l) : 
lMiller&ScaloHl979l) .l Kuruczl dl992T) atmosphere models , for pop- 
ulation synthesis, and Padova tracks dBertelli et alii 19941) are also 
considered. Moreover, SSPs of GRASIL cover a spectral range from 
91 A up to 1200/im. Further details can be found in Laura Silva PhD 
thesis 



2.3 Maraston models (M05) 

The Maraston EPS modelsQ (he reafter M05) are being developed 
by Cl audia Maraston since 1998 dMarastonl 19981) with an update in 
2005 dMarastonl2005l) . They are based on the fuel consumption the- 
orem and include a proper treatment of the TP-AGB phase. Accord- 
ing to these models, the effects of TP-AGB stars in the NIR spectra 
are unavoidable. The M05 model s, by including empiric al spectra 
of oxygen- and carbon-rich stars dLancon & Woodll200d) , can pre- 
dict t he presence of NIR absorption features such as the 1.1 /im CN 
band dRiffel et alj2007il2008l,l2009l) . whose detection can be taken 
as an unambiguous evidence of a young to intermediate age SP. The 
models have been used by our team to study stellar populations 
in active galactic nuclei and starburst galaxies (Riff el et al]|2007l . 
1 200 8ll 20091. [2010l. for CN, see also Ramos-Almeida et al. 2008 and 
Dottori et al. 2005) as well as in the age dating of massive galax- 
ies at high redshjft dMaraston etaill 20061: Ivan der Wei et alj|2006l : 
iRodighiero et alj|2007l ; lcimatti et alj2008^ 

M05 models span a range of 6 different metallicities (^j ^ 
^3.5) with ages distributed from 1 Myr to 15 Gyr according 
to a grid of 67 models (Note that the full age grid is n ot avail- 
able fo r all metallicities, M 05). The IMFs considered are: ISalpetej 
dl955l) and lKroupdd200ll) . The stellar spectra were also taken from 
the BaSel library. The spectral range is from 9lA to 160/im, with 
a spectral a resoluti on of 5-10A up to the o ptical region, and 20- 
100A in the NIR. In lMarastorJ dl998ll2005h models, the TP-AGB 
contributes with 40% to the bolometric flux, but rising to 80% when 
only the Jf-band is considered. 



2.2 GALEV 

The GALaxy EVolution - GALE\0 - evolutionary synthesis models 
describe the evolution of stellar populations in general, from star 



2.4 GALAXEV 

GALAXE\J3 is a widely used library of evolutionary stellar popu- 
lation synthesis models. It is computed with the isochrones syn- 
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thesis code of lBruzual & Charlotl d2003l also known as BC03) The 
spectral coverage of this library is from 91 A up to 160/im, with 
a resolution of 3 A between 3200 and 9500 A, and a lower res- 
olution elsewhere. Ages range from 1 x 10 5 up to 2 x 10 10 yr, 
for a wide range of metallicities (»L < -£- < 2.5). These 
mode l s use the STELIB/B aSeL libraries (see lLeieune et alJI 19971 
Il999l : IWestera et alj|2002l. and referenc es therein) as well as the 
STELIB/Pickles libraries dPicklesll 19981). GALAXEV allows the use 



of two IMFs jChabrierll2003l; ISalpeteJj 19551) and 3 stellar evolu- 
tion t r acks: Geneva i Schaller et al.ll 1992^). Padova 94 (lAlongi et al. 
19931 : iBressan et al.l 1 19931; iFagotto et all Il994tj|bl iGirardi et al 



1996) and Padova 00 iGirardi et al. 2000]). These models do not in- 



clude the TP-AGB phase. 



the 5Myr population at 3800A represents only 5% at 1.223/im). 
However, the intermediate to old components (t > 500 Myrs) can 
be directly compared between different wavelengths. This does not 
occur when dealing with mass fractions, since the age derived by 
the mass fraction is a more physical parameter, but has a much 
less direct relation with the observables, depending strongly on the 
M/L ratio, which is not constant. Thus, the light fraction can be 
taken as an direct observable parameter and use Eqs.[TJand[2] to 
propagate the results over other spectral regions (see AppendixlAl. 

However, the stellar population of galaxies is not so simple as 
two components. Therefore, we have divided our SSPs into 3 popu- 
lation vectors a-j^O.005 Gyr; Xi=0.025...1.0 Gyr and x =13 Gyr. To 
investigate the effect of adding one more component to the above 
exercise, we combined three population vectors according to: 



3 METHODOLOGY AND RESULTS 

First, we investigate the dependence of the stellar population 
components on the normalisation point, from the NUV to the 
NIR. To do this, we select spectral regions free from emis- 
sion/absorption lines to be used as normalisation points, Fa. They 
are: 3800A 4020A 4570A 5300A 5545A 5650A 5800A 5870A 
6170A 6620A 8100A 8815 A 9940A 1 058/xm 1.223»m 1 520^m 
1.701jnm 2.092/jm 2A9ym llBicall 19881. iRickes et alj|2008l. 120091 : 
ICid Fernandes et alj2004l : ISaraiva et al]200ll ; lRiffel et alj2008ir 
In addition, we select SSPs with ages: 0.005, 0.025, 0.050, 0.1, 
0.2, 0.3,0.4 0.5,0.6, 0.7, 0.8, 0.9, 1.0, 2.0, 5.0 and 13 Gyr as repre- 
sentative of the stellar populations observed in galaxies. As a first 
exercise we combine two components: 13 Gry (the old population) 
and one of the other SSPs representing the "young" population. 
The combination was made by summing up, along the whole spec- 
tral range (~3500A to 2.5^im), increasing fractions of the "young" 
component from 1 to 100%, according to 



F=(l-f)F y + fF ; 



(1) 



where / is the fractional flux, which we vary in steps of 0.01, F y is 
the flux of spectrum of the "young" component, for each A between 
~3500A to 2.5/im, normalised to unity at 5870A for the optical, or 
at 1.223/im for the NIR, and F is the flux of the 13 Gyr spectrum 
also normalised at the same points. Note that the normalisation of 
the spectra is done by dividing their fluxes by that of the normali- 
sation point (5870A or 1.223/^m), after the computations are done. 
In addition, by normalising the spectra we are dealing with light- 
fractions, which are directly related to the observations. F is the 
resulting flux of the combined SSPs (young + old) at a specific 
A. The averaged stellar population components spread over all F\ 
were derived by: 



(2) 



Eq|2]represents the young and old light fractions for different 
Fx, this is what we call fraction at A. The result of such process 
is summarised in Figs. [T] and [2] In these figures we show the sum 
of a "young" component with a 13 Gyr SSP (the old component) 
in steps of 10%. In practice we start with 5% of the flux of the 
"young" SSP + 95% of the flux of the 13 Gyr SSP, and finish with 
95% of the young + 5% of the old component. This procedure is 
done over all As between the NUV and NIR. It is worth mentioning 
that Figs. [T] and [2] show M05 EPS models as reference (Sec.[4]l. 
These figures suggest that one cannot directly compare the light 
fractions of the young component (t < 400 Myrs) derived in the 
optical with those obtained in the NIR, and vice versa (i.e. 77% of 



F ' = Too (7Fy + 5Fi + 77Fo); 

with 7 + 5 + 77 = 100 



(3) 



where, 7, S and 77 are the fractional fluxes from 0% to 100%, which 
are varied in steps of 1%; F y , Fi and F are the normalised fluxes 
(at A=5870A or A 1.223^im) of the x y , x\ and x population vec- 
tors, respectively. F' is the resulting flux of the combined SSPs 
(young + intermediate + old) at a specific A. 

The averaged stellar population components distributed over 
all As were derived similarly to Eq.[2] 



Jy jp/ ; J 1 rp, Jo 



(4) 



We show the result of such combination in Figs. [3] to [5] It is 
clear that even with three population vectors, the results derived 
in one wavelength are not the same as in other As. In addition, 
Fig.|5]reinforces the fact that in the case of only intermediate to old 
components, the population fractions derived in one wavelength are 
nearly the same as in other As (i.e. the values derived in the optical 
can be used in the NIR). 



4 DISCUSSION 

The determination of mean ages and metallicities of galax- 
ies is important for models of galaxy formation and evolution. 
The techniques for dating unresolved stellar populations in local 



sured in the NUV and optical (e.g. 


Bidll988l: IWorthevlll994 


Bica et al. 19961: Bonatto et al. 19961. 


19981; Traaer et all 2000dlb 


Maraston & Thomas 2000; Thomas et al. 2005; Rickes et al. 2008 



and references therein). However, even small mass fractions of 
young stars added to an old population can affect the NUV and 
optical age determinations significantly, making the galaxy appear 
young, whic h leads to a further degeneracy between ma ss frac- 
tion and age ^Thomas & Daviesl2006l : ls"erra & Tragej2007l) . These 
problems may be solved for intermediate-age stellar populations 
looking in NIR wavelengths. In addition, the detecti on of young 
stella r populations in the NIR requires a hard work jRiffel et al .1 
l20ld) . Nevertheless, it is clear from Fig. [2] that even a small frac- 
tion of a 5 Myr population (~ 5%) detected in the NIR may be 
responsible for almost all the light observed in the NUV (~ 70%). 

Clearly, synthesis results should not be directly propagated 
from the NIR to the NUV/Optical, or vice versa. Instead, Eqs. [2] 
to |4]should be used for this purpose. To help with such a compar- 
ison we have created an on-line form, the Panchromatic Averaged 
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Figure 1. Light-fraction at A for the stellar population components. The curves result from Eq[2]with increasing fractions of 10% of the "young" component 
from bottom to top. In each panel we start with 5% of the flux of the "young" SSP + 95% of the 13 Gyr SSP, and finish with 95% of the young + 5% of the 
old component. The ages of the "young" component are on the labels. The normalisation point for the combination was 5870A. We use the M05 models as 
reference. See text for more details. 



Stellar Population: PaASiS and make available for download the 
tables with the results of the above equations (see Appendix [At. 
Another important ingredient in stellar population fitting is the 



available at: http://www.if.ufrgs.br/~riffel/software.html 



metallicities used. As shown by IChen et al the results of 

the fitting have a weaker dependence on metallicity than age. The 
question which arises here is, does metallicity affect the propaga- 
tion of the averaged stellar populations? We investigate this effect 
with M05 SSPs with 3 different metallicities (i^ZO; ZQ; and 
2 ZQ) and the same age grid as in Fig. [T] The results are shown 
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Figure 2. Same as Fig.[T]for the normalisation point at 1.223/^m. 



in Figs. [6] and [7] Clearly the propagation of the contributions has 
a negligible dependence on metallicity. Thus, one can use the con- 
dense d population vectors proposed bvlcid Fernandes et al, U2004 
to propagate the fitting results over all As. 

All the tests described above were made using the M05 mod- 
els, but, as stated in Sec. [2] there are more EPS models available in 
the literature covering simultaneously the spectral region between 
~3500A and 2.5/im. Thus, it is necessary to test if the selection 



of EPS models will produce different results in the propagation of 
the synthesis results over different As. In Fig. [8] we compare the 
different models among each other. It is clear that in the case of 
the optical normalisation point (5870A), the four models produce 
very similar results in the interval between 3800 A and 9000 A, but 
a discrepancy between GALEV/M05 and GRASIL/BC03 models is 
observed in the NIR. Such a discrepancy is due to the well known 
fact that GALEV and M05 models do include stars in the TP-AGB 
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Figure 3. Light-fraction at A for the stellar population components. The curves result from Eq|4]with increasing fractions from bottom to top in each panel. 
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Figure 4. Same as Fig.[5]for a different set of ages. 
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Figure 5. Same as Fig.|5]for a different set of ages. 



phase (see M05, for example), which is more sensitive to the NIR 
than the optical, i.e TP-AGB stars account for 25 to 40% of the 
bolometric light of an SSP, and for 40 to 60% of the light emitted in 
the K-band (see lSchulz et al.ll2002l ; lMarastonll2005l and references 
therein). However, there is a difference between GALEV and M05 
models, enhanced in the 100 Myr population. There are two pos- 
sible explanations for this discrepancy: one is associated with the 
different onset age of the TP-AGB o n the models. A high TP-AGB 
contribution at 100 Myr, as applied bv lSchulz et al. I l2002h , which is 
excessively high when compared t o young Large Mag ellanic Cloud 
globular clusters jMarastoJl99llMarigo et al.l2008l) . The other is 
associated with the w a y in which the TP-AGB treatment is made 
dMaraston etal] 12004 lBruzualll2007l) . GALEV includes TP-AGB 
by means of isoch r ones (Padova94 + i mproved TP-AGB models 
iBertelli et alj| 19941 ; iGirardi et alJl200d : lMarigo et alj|2008l) , while 
M05 is based on a different approa ch, the fuel consumption theo- 
rem. According to iMarastonl ( 120051) ■ the stellar luminosity during 
the evolutionary phases that follow or suffer from mass loss can- 
not be predicted by stellar tracks, because there is no theory linking 
mass-loss rates to the basic stellar parameters, such as luminosity. 

Note that the trend observed in Fig. [8] can also be associated 
with the fact that in M05 models, the TP-AGB contributes with 
40% to the bolometric flux, and 80% to the K-band. This is higher 
than the "simp ler" calculati o ns (m ade by the Padova group at the 
time) used by ISchulz et aljj 20021), in which only some thermal 
pulses have been included I Girardi et alj|2000l) . 

In addition, the tests were performed with the lSalpeteil Jl955h 
IMF. To see the effect of the IMF on the averaged stellar popula- 
tions over all As, we repeat the same exercises for different IMFs 
and summarise the results in Fig. [9] This figure suggests that the 



IMF does not play an important role in the synthesis propagation to 
other spectral regions. 



4.1 Testing PaASP 

In order to test if we are able to predict the averaged stellar popu- 
lations over all lambdas, extensive simulations were performed to 
evaluate the PaASP's ability to recover inp ut paramet e rs. In short, 
we build artificial galaxy spectra by mixing IMarastonl J2005f) solar 
metallicity SSP models of 5 Myr, 700 Myr and 13 Gyr in different 
proportions, and perform stellar population synthesis with normali- 
sation points at the optical (5870A) and NIR (12230A). To this pur- 
pose we use the STARLIGHT cod e optimised to the optical region 
jCid Fernandes et al.|[2004l 120051 : iMateus et alj|2006t lAsari etal] 
l2007l : ICid Fernandes et alj2008t) . and wh ich also produces reliab le 
results when applied to the NIR spectral dRiffel et al .11 20091 120 1 Oh ■ 
In summary, STARLIGHT fits an observed spectum Ox with 
a combination, in different pr oportions, of N* SSPs. Ba sically, it 
solves the following equation idd Fernandes et al]|2005l) : 



' -2V* 

i=i 



<g> G(u*,er») 



(5) 



where M\ is a model spectrum, bj,\r\ is the reddened spectrum 
of the jth SSP normalised at A ; r x = 10 _0 ' 4(Aa " Aao) is the red- 
dening term; Mao is the synthetic flux at the normalisation wave- 
length; x is the population vector; ® denotes the convolution oper- 
ator, and G(v*, <t*) is the Gaussian distribution used to model the 
line-of-sight stellar motions, which is centred at velocity v* with 
dispersion a, . The final fit is carried out with a simulated anneal- 
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Figure 6. Metallicity effect for the normalisation point at 5870A. Boxes represent solar metallicity, circles jgZQ and triangles are 2 ZQ, We use the M05 
models as reference. 



ing plus Metropolis scheme, which searches for the minimum of 
the equation: 

X 2 = Yltt°* ~ M >>a] 2 , (6) 

A 

where emission lines and spurious features are masked out 
by fixing W\=0. For a detailed description of STARLIGHT see 
ICid Fernandes et al.|j2004l . |2005h . 



As base set we take lMarastonl d2005l) SSPs covering 14 ages, 
t= 0.001, 0.005, 0.01, 0.03, 0.05, 0.1, 0.2, 0.5, 0.7, 1, 2, 5, 9, 13 Gyr, 
and 4 metallicities, namely: Z= 0.02 Zq, 0.5 Zq, 1 Zq and 2 Z@, 
summing up 56 elements. 

To compare STARLIGHT with PaASP predictions, we use the 
condensed population vector, which is obtained by binning the syn- 
thesis results into young, Xy (t ^ 5 x 10 7 yr); intermediate-age, xj 
(1 x 10 8 < t < 2 x 10 9 yr) and old, xo (t > 2 x 10 9 yr) compo- 
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Figure 7. Same as Fig.|6]for the normalisation point at 1.223/^m. 

nents jRiffel et alj2009tlCid Femandes et alj2004h . These compo- all points fall in a region less than 5% from the identity line). The 
nents were then taken to represent the 5 Myr, 700 Myr and 13 Gyr confidence level drops to ~85% in the NIR. 
old populations. These vectors are compared with the PaASP pre- 
dictions in Fig. 1101 

5 FINAL REMARKS 



Clearly, our predictions are consistent with the stellar popu- 
lation synthesis, especially at the optical region, where the confi- 
dence level between predictions and synthesis is ~95% (i.e. almost 



We study the panchromatic stellar population components over the 
3500A to 2.5/im spectral region. In particular, we analyse how 
the spectral fitting of galaxies based on light-fractions derived 
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Figure 8. Comparison of different EPS models. Normalised at 5870A (left) and 1 .223/im(right). Note that we use 100 Myr + 13 Gyr SSPs. 



in a given spectral range can be propagated over all As. De- 
pendencies on EPS models, age, metallicity, and stellar evolu- 
tion tracks of four widely used EPS models (grasil, galev, 
Maraston and galaxev) were taken into account. Our main re- 
sults are: 

• The young (t < 400 Myr) stellar population fractions derived 
in the optical cannot be directly compared to those derived in the 



NIR, and vice versa. For example, a contribution of ~ 80% of a 
5 Myr population at 3800A translates into only 5% at 1.223/im. 

• The intermediate to old age components (t > 500 Myr) can be 
directly compared from the NUV up to the NIR. 

• The metallicity dependence on the propagation of the stellar 
population fractions derived from NUV — NIR is negligible. 

• Different EPS models produce similar results in the propa- 
gation of the synthesis results over the interval between 3800 A 
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Figure 9. Effect of different IMFs. At left the normalisation point is 5870A and at right 1.223/im. Salpeter (SP) and Kroupa (KP) where taken from M05 
models; Kennicutt (Ken), Scalo (Sea) and Miller & Scalo (MSc) IMFs were taken from GRASIL EPS and Chabrier (Chab) are from GALAXEV. 
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Figure 10. PaASP tests. Stars, boxes and circles represent the 5 Myr, 700 Myr and 13 Gyr SSPs, respectively. Identity is shown by the solid line and dotted 
lines indicate the 90% confidence level. 
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and 9000 A. However, a discrepancy between GALEV/M05 and 
GRASIL/BC03 models occurs in the NIR. Such a discrepancy may 
be due to the fact that GALEV and M05 models do include stars in 
the TP-AGB phase. 

• There is a difference between GALEV and M05 models, en- 
hanced in the 100 Myr population. Such a discrepancy may be 
associated with the way in which th e TP-AGB treatment is made 
dMaraston etaifeoOrJ: lBruzualll2007l) . 

• We test the effect of 6 different IMFs, and we conclude that 
the IMF is not important in the propagation of the synthesis results. 

• Extensive simulations were performed to evaluate PaASP's 
ability to recover input parameters. Our predictions are consistent 
with the stellar population synthesis with a confidence level be- 
tween the predictions and synthesis of ~95% in the optical and 
-85% in the NIR. 

In summary, spectral fitting results should not be di- 
rectly propagated from the NIR to the NUV/Optical, or vice 
versa. Instead, Eqs. [2] to [4] should be used for this pur- 
pose. However, since this is hard to do when dealing with 
large samples of objects, we have created an on-line form, the 
Panchromatic Averaged Stellar Population - PaASP, available at: 
http://www.if.ufrgs.br/~riffel/software.html We also make avail- 
able for download the tables with the results of the above equations 
for a wide range of ages (see Appendix [At . 
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APPENDIX A: HOW TO USE THE PaASP FORM 

In this appendix we provide details how to use the PaASP form. 
PaASP is available at http://www.if.ufrgs.br/-~riffel/software.html 

(i) Bin your synthesis results (Xj) into 3 population vectors: 
young, intermediate, and old ages. Fo r more information on the 
defini tion of the populations v ectors see lCid Fernandes et al. I d2004l 
l2005h and lRiffel et all J2009h . 

(ii) Select the representative age of each vector, for example: 
5 Myr for young, 200 Myr for intermediate, and 13 Gyr for old. Tip: 
take as the representative population the ages of Xj with the largest 
contribution in each bin. 

(iii) Put the percentage contribution (only integers and sum 
equal to 100% are allowed) of each age in the form, select the nor- 
malisation point and submit it. 

(iv) You will be redirected to the results query page. There your 
inputs are marked in red and the propagated results over all As (see 
text) are shown. You can also download the file with a table with 
propagation results for all the results of Eq.[4] in fractions of 1% 
for your query. 

(v) The 3 first columns of the table are the input fractions for a 
given normalisation point and for 3 ages, which are specified in the 
header (first and second lines). The next columns are the propaga- 
tion of the results for all As. Note that the results are given in 3 lines 
blocks: young, intermediate and old fractions, respectively. 

PaASP is freely distributed and is supported by Brazilian fund- 
ing agencies CAPES and CNPq. An acknowledgement for the use 
would be appreciated. 



